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Abstract 

We report on the computation of the angle-between-leptons distribution in Drell-Yan processes. 
More precisely we study the recently introduced variable 4>* , which provides us with a more accurate 
probe of the low Qt domain of Z boson production at hadron colliders. Our theoretical prediction 
is obtained by matching a next-to-next-to leading logarithmic (NNLL) resummation to a fixed order 
calculation at next-to-leading order (NLO). We find that the result significantly differs from a pure fixed 
order calculation in a wide range of the observable we are interested in, clearly indicating the need for 
resummation. We also perform a first comparison to the measurement by the D0 collaboration, finding 
good agreement between theory and experiment. 
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1 Introduction 



The production of a lepton pair in hadron-hadron collisions is one of the most studied processes in 
particle phenomenology, with the original paper appearing more than forty years ago pQ. Since then 
a huge theoretical effort has gone into improving the accuracy of the predictions. For instance, QCD 
corrections are known to 0{a 2 s ) In particular, the transverse momentum distribution of the lepton 
pair, or equivalently of the gauge boson decaying into it, is of great interest. It is sensitive to multi- 
gluon emission from the initial state partons. This is a classical example of a multi-scale problem and 
the correct treatment of these effects goes beyond fixed order perturbation theory. Let us introduce 
the invariant mass of the leptons M, that will be chosen around the Z mass, and the let us call Qt 
the magnitude of the Z transverse momentum. We have to consider three different regimes. When 
Qt ~ M we expect fixed order perturbation theory to work and programs like MCFM [3] will give a 
good description of the process. In the region Aqcd "C Qt « M, we can still rely on perturbation 
theory but large logarithms of the ratio Qt/M may spoil the convergence of the perturbative expansion 
and must be resummed to all orders. Finally, in the region Qt ~ Aqcd we expect non-perturbative 
effects to play a significant role. Therefore, it is important to compute a solid perturbative prediction, 
so that we can compare it to precise data coming from the experiments and be able to pin down non- 
perturbative contributions, related to the intrinsic transverse momentum of the initial-state quarks. 
An accurate theoretical description of the transverse momentum of the Z boson is also relevant for the 
extraction of the W mass. 

The resummation of large logarithms in the Qt spectrum has been studied by several groups for 
many years and it is currently known to NNLL accuracy (see Ref |2] and references therein). Accurate 
theoretical predictions have been compared to data coming from the D0 and CDF experiments at the 
Tevatron but no clear conclusions have been drawn to date about the relevance of non-perturbative 
effects. One of the limiting factor has been the experimental resolution which affects the measurement 
of transverse momenta. For this reason novel variables have been introduced in [5j[6] and recently 
measured by the D0 collaboration [7j. These variables, labelled the and eft*, both crucially depend 
on the azimuthal angle A(ft between the final state leptons, at low Qt- The experimental resolution for 
aT and <ft* is significantly better than the one for Qt [6]. 

The D0 collaboration compared the result of their measurement to the theoretical prediction of 
the program RESBOS [8 J . They found overall agreement with some discrepancies in the large rapidity 
region. In particular, the data disfavour current non-perturbative models, such as small- £ broadening [9] 
. Therefore we need an accurate theoretical prediction for these new variables, along the line of those 
for Qt resummation, to be able to assess the importance of non-perturbative effects. 

2 The (J)* distribution 

In Ref. [10] we have computed a theoretical prediction for the (ft* distribution by matching a resummed 
NNLL calculation, which captures the dominant behaviour at small (ft* to a NLO one obtained from 



1 



the program MCFM: 

da \ __ (do_\ f (—} - ( — 

^0* J matched V ^0* / rcsummcd V / NLO V ^0* / expanded 

The resummed distribution has the following form 

^-(ct>*,M,cosO*,y) = ^ [°° dbM cos {bM<j>*)e- R & M >^>^ (2) 
d(f>* sN c Jq 

x S (x 1 , x 2 , cos 9*, b, M, fJ,Q,fJ,R, /xp) , 

where xi^ = "^= e±?/ an d 6 = ^p-- The above result is yet to be integrated over the dilepton invariant 
mass M, the scattering angle 0* and rapidity of the dilepton system (or equivalently the Z boson rapid- 
ity) y. Note the dependence on three arbitrary scales: renormalisation, factorisation and resummation, 
which can be varied in order to estimate the theoretical uncertainty. At the moment we set them all 
equal to each other and to the dilepton invariant mass M. 

The dependence upon the large logarithms we wish to resum is encoded in the radiator: 

R (bM) = LgW (a s L) + g® (a s L) + %< 3 ) (a s L) , (3) 

7T 

where L = ln(6 2 M 2 ) and a s = a s (M). The functions g^> are the same as in [3] but in our case 
the radiator does not contain any term involving the DGLAP anomalous dimensions or the coefficient 
functions. These iV-dependent contributions have been used to evolve the parton distribution functions 
or have been taken care of by evaluating the running coupling in front of the coefficient functions 
at the appropriate scale 1/6. The coefficient appearing in the function g% has been recently 
determined [11] . However, we currently include in g^ only the terms which are relevant at O (a 2 ), and 
we call our partial NNLL resummation NNLL*. 

Before presenting results for the matched distribution Eq. (TT]) we must check that the expansion of 
our resummation to O(a^) agrees with the NLO calculation. We remind the reader that the relation 
between the resummation for Qt and ay was worked out in [12] and we also have that eft* ~ a^jM 
as small Qt- In order to check our understanding of the relation between the different observables we 
compute 
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AD{e) = V dln7 ^ ^ ~ - (Gr/2) I Qt/2 =J • (4) 

Subtracting AD from the corresponding fixed-order differential distribution D (<f>*) — D (Qt/2), com- 
puted with MCFM at NLO, we find a result that tends to zero, as shown in Fig. [H on the left. Thus, 
we have complete control over the divergent pieces at order c? s and we can adopt a particularly simple 
matching formula like the one in Eq. ([1]) , where one adds the resummed differential distribution to the 
NLO result and subtracts the expansion of the resummation to order a 2 . 

The result for the matched differential (jf distribution is plotted in Fig. [1] on the right, together 
with the pure fixed order O (a 2 ) calculation obtained from MCFM. The curves are obtained taking into 
account the D0 cuts for the muons, i.e. 70 < M < 110 GeV, pt > 15 GeV and rj < 2, and integrating 
over the rapidity of the Z boson. We notice that the NLO calculation diverges in the region of small 
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L=ln(e) 4>* 



Figure 1: On the left: the difference between the NLO differential distributions for <jf and Qt/2 from 
MCFM after removal of logarithmic terms from the resummation. On the right: the differential <j? 
distribution computed at NLO with MCFM (dashed blue line) and our final result (solid red) obtained 
by matching the NNLL* resummation to the fixed order calculation at 0(a^), from MCFM. The 
distributions are normalised to NLO cross section. 



<j>*, while the resummed and matched results tend to a constant, without forming any peak. While low 
values of <p* or Qt can be obtained via Sudakov suppression or kinematical cancellation, in the present 
case of <p* the kinematical cancellation starts to dominate prior to the formation of the Sudakov peak, 
in contrast to the Qt case. 

We are now in the position to perform a first comparison to the data obtained by the D0 collabo- 
ration [7j. For the moment we restrict ourselves to the muon case in two different rapidity bins: \y\ < 1 
in Fig. [2j on the left and 1 < \y\ < 2, on the right. We can clearly see that in both cases the matched 
result gives a good description of the data, with a maximum discrepancy of 7/8 % in the last (f>* point. 
This is a very good result because our curves are obtained from a pure perturbative QCD calculation, 
with no other ingredients. The insets show that the matching has beneficial effects also at large (ft*. 

The plots in Fig. [2] are very encouraging but still not complete. First of all we have to include 
the complete expression for the function g%. Moreover, before making any statement about the size on 
non-pertubative effects, we have to study the theoretical uncertainty. More precisely we can make an 
estimate of missing higher orders in the fixed-order part by varying renormalisation and factorisation 
scales and hf- We can also estimate missing higher logarithmic orders by varying the resummation 
scale hq. Moreover we need to investigate different sets of parton distribution functions and different 
procedures to regularise the 6-integral. Having done that, we will be able to add an uncertainty band 
to our curves, properly compare to the D0 data and pin down the non-perturbative contribution. 

From a more theoretical viewpoint we note the resummation presented here is closely related to the 
one for Acp between jets |13j . which therefore signifies an extension of Qt resummation to processes 
with colour in the final state. 
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Figure 2: Comparison of our resummed and matched theoretical prediction to the muon data collected 
by the D0 collaboration. The plot on the left corresponds to the central rapidity bin \y\ < 1, while the 
one on the right is for 1 < \y\ < 2. The curves have been obtained without non-perturbative effects. 



3 Conclusions 

We have reported on phenomenological work on accurate predictions for novel variables recently mea- 
sured by the D0 collaboration to probe the low Qt region of the Z boson spectrum. Our result resums 
large logarithms at NNLL* accuracy and it is matched to a NLO calculation. The first comparisons 
to data are encouraging, although before drawing any conclusions about the size of non-perturbative 
effects, we must estimate the theoretical uncertainty. 
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